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Abstract 

We propose an expression for the current form of the lowering operator of 
the sl2 loop algebra symmetry of the six vertex model (XXZ spin chain) at 
roots of unity. This operator has poles which correspond to the evaluation 
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From these polynomials we find that the Bethe roots satisfy sum rules for 
each value of S z . 
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I. INTRODUCTION 



Recently it has been demonstrated@-i that the six vertex model with periodic boundary 
conditions specified by a 2 L x 2 L transfer matrix T(v; 7) 

T {jl}t{jlf} {v;i) = £ ^ 1 y i (A 1 ,A 2 )W M2y2 (A 2 ,A 3 )---^yjA £ ,A 1 ) (1.1) 



A i= ±l 



with 



W^{\i, Ai+i)= -(a + 6)^,^5^,^! + «(o - b )^fl,x i+1 



where cx J are the Pauli matrices and 



+c(a + ,cr7 , +o- _ ,o-t i ). (1.2) 



a = isinh -(1; — 27), 6 = — i sinh —(v + ij), c=— ismhij (1-3) 
has an s/ 2 loop algebra symmetry when 7 satisfies the "root of unity condition" 

with m and iV relatively prime. The operator 

s z = \i:^ (1.5) 

commutes with the transfer matrix T(t>;7) and in the sector S* 2 = 0(mod N) the sZ 2 loop 
algebra symmetry was proven in ref.0 by demonstrating that the operators with q = e^ 71 " -70 ) 

S ±(N) = £ q (N/2)«* q . . . 3 g (W g ^ ^ q W-W)o» g . . . 

l<n <•••<?><£ 

g ? ((iV-2)/2) CT * ^ a ± ^ g ((iV-4)/2)^ ^ . . . g a ± g g -(iY/2 K g . . . g, g ~(W (l g) 

and 

= ]T <T W2)ff * ® • • • ® ^ W2)CTZ ® ® g- ((Ar ~ 2)/: V ® ■ • ■ 
i<ii<-<ijv<i 

g ff ± g g . . . g, o-± g, ? (W g . . . g, ? (W (1.7) 

commute (anticommute)with the transfer matrix if iV — m is even (odd) and satisfy the 
defining relations for the Chevalley generators of the sZ 2 loop algebra. In the sectors where 
S z = r 7^ (mod N) with 1 < r < N — 1 operators such as (T + ) r (S , ~) r S'~( iV ) commute 
(anticommute) with the transfer matrix and the sl 2 loop algebra symmetry was inferred 
from computer computations. 

The eigenvalues of the transfer matrix are polynomials in the variable e v . In the limit 
v — ► ±^7 the transfer matrix reduces to 
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T(v;y) — > (sinhz7) L n± 



/ - —r 1 - H xxz + —-L + 0((v =F «7 2 
smnry V 2 / 



;i.8) 



where Il± is the left (right) shift operator whose eigenvalues are e TlP where P is the mo- 
mentum of the state and 



with 

A = -cos7 (1.10) 

is the Hamiltonian of the XXZ spin chain. Therefore every degeneracy of the transfer 
matrix (polynomial in e v ) eigenvalues implies a corresponding degeneracy of the (numerical) 
Hamiltonian eigenvalues. However, the converse need not be true and the Hamiltonian ( jl.9| ) 
is known to have degeneracies such as that caused by invariance under parity which the 
transfer matrix does not have. 

The existence of this symmetry means that the space of eigenstates states of the transfer 
matrix can be decomposed into finite dimensional representations of the si 2 loop algebra. 
Each finite dimensional representation has a vector \Q > with the property that 

s + w\n>=T + ^\n>=o (1.11) 



and we define S* by 



s z \n >= s: a jn > . (1.12) 

In analogy with the finite Lie algebra sl% we call this (following ref.i) the highest weight 
vector of the representation. 

A fundamental property of all affine Lie algebras^ is that they may be defined either by 
a Chevalley basis or by a mode basis. In the mode basis the elements of the sl% loop algebra 
are e(n), f(n) and h(n), where n is an integer, which satisfy the commutation relations 

[e(m), f(n)}= him, + n) 
[e(m), h(n)]= —2e(m + n) 

[f(m),h{n)]=2f{m + n) (1.13) 
and there is the relation to the Chevalley basis of 

e (0)=T-W e(-l) = S-W 

/(0) = T +(7V) , f(l) = S +{N) . (1.14) 

The theory of finite dimensional representations of affine Lie algebras has been extensively 
studied in!-! in terms of this mode basis where it is shown that all irreducible finite di- 
mensional representations are tensor products evaluation representations. These evaluation 
representations are specified by vectors \aj,m,j > where for all integer n (positive, negative 
or zero) 
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e(n)\aj,rrij > = Oje mj \aj,mj > 
/(n)|aj-,mj > = " •./',,., a r m s > 
h(n)\aj,rrij > = a^h mj \aj,rrij > 



;i.i5) 



where are called evaluation parameters and e m . , f m . , h m , are a spin m,j representation of 
5/2- The evaluation parameters are shown to be the roots of the Drinfeld polynomial Pn(z) 
defined as 

Pn{z) =\{{z - a 3 )^ (1.16) 
3 

where the aj are distinct and 

^(^ElirK (1.17) 
r>0 

and for each highest weight state \Q > we compute \i r from the eigenvalue equation 

rp+(N)r g-(N)r 



\tt>=fi r \tt>. (1.18) 



If the roots of the Drinfeld polynomial ( |1.17| ) are all distinct only spin 1/2 representations 



occur. In this case if we call dn the degree of the Drinfeld polynomial then the number 
of evaluation parameters is d n and the number of states with spin S z = S^ax — ^ ^ s ^ ne 
binomial coefficient 




= number of degenerate eigenvalues (1-19) 

In the mode basis we define the generating function (currents) for sufficiently large z by 

00 

E-(z) = Y J <n)z~ n . (1.20) 

n=0 

From the definition of evaluation representation this current will have poles when z = aj. In 
order to show that only spin 1/2 representations occur it is sufficient to demonstrate that 
there are only single poles in the current E(z). 

To compute the Drinfeld polynomial from the definition Ql.17 ) we need to be in possession 



of the highest weight vectors. For a finite Lie algebra these highest weight vectors cannot be 
computed from symmetry considerations alone. However for affine Lie algebras if the sets of 
evaluation parameters are all distinct then the algebra will be powerful enough to determine 
the highest weight vectors for all multiplets except singlets. 

For the six vertex model, however, we may make use of the fact noted in ref.sBJ! that 
the highest weight vectors of the si 2 loop algebra are identical with the eigenvectors which 
are obtained from the solution of Bethe's equation for j = 1, • • • , ^ — \S Z \ 




Vi + 2ry) 

9 . 1.21 



4 



In fact there are two distinct ways of specifying the eigenvectors in terms of the solutions 
(Bethe roots) of the equation Q1.21 ); the coordinate Bethe ansatz§-lll which uses explicit 
forms of the wavefunctions in coordinate space and the algebraic Bethe ansatJU. which 
produces the wave functions by exhibiting an operator which creates the wave functions by 
acting on the state of all spins up. 

In this paper we use the algebraic Bethe ansatz. In this method of solution of the six 
vertex model we introduce the (2 x 2) monodromy matrix (of operators) 

Mx,y(v; 1 )= £ W Wrf4 (A, A 2 )W^ 2 (A 2) As) ■ • ■ W^Xl, A') (1.22) 

A i = ±l 

for which we will use the conventional notation@-0 

' A{v;i) B(i 
C(v;y) D(v;j) 



The transfer matrix (|1.1|) is 

T{v- 1 )=A{v- 1 ) + D{v- 1 ). (1.24) 

and here and in the rest of the paper we make (for convenience) the restriction that L be 
even and thus the eigenvalues of S z are integer. 

Fundamental to the algebraic Bethe's ansatz are the "commutation" relations of the 
operators (2 L x 2 L matrices) v4(t>;7), B(v;j), C(v]j), D(v;j) of which for our purposes 
here we will need the following three 

[B{v),B(v')}=0 (1.25) 

A{v)B{v')= f(v - v')B(v')A(v) + g{y' - v)B(v)A(v') (1.26) 

D(v)B(v')= f(v' - v)B(v')D(v) + g(v - v')B{v)D{v') (1.27) 

where 

. , sinh h(v + 2z7) 

f( v )= 2 ,\ ; - u 1.28 

J v ; sinh(w/2) v ; 

. , sinhi7 

9(v = - . . 1.29 

smh(f /2) 

and note the periodicity properties 

B(v -2m) = -B(v), C{v-2m) = -C(v) } A(v - 2m) = A(v), D(v - 2m) = D(v). 

(1.30) 

The eigenvectors of the six vertex model for all states where the solutions to Bethe's 
equation (|1.21| ) do not lead to undefined factors of 0/0 are given in this notation by 

n 

\{v k }>= Y[B{v k - n )\0> (1.31) 

k=l 
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where |0 > is the unique state with all spins up with 



and 



S 2 |0>=^|0> (1.32) 



S z \{v k } >= {\ - n)|K} > • (1.33) 



The corresponding eigenvalue of the transfer matrix is 

\ n 1 71 

t(v) = sinh L -{v - ij) J] f(v - Vj) + sinh L -{v + i 7 ) [] /(% " (1-34) 

^ i=i 1 i=i 



If 7 is a "generic" value then the eigenstates ( 1-31 ) are complete0-^l and there are no 



factors of 0/0 in Bethe's equation ( J1.21 ) . However, when the root of unity condition ( |1.4| ) 
holds it was seen on refs.ELi that factors of 0/0 do occur in ( [L21| ) because of solutions which 
approach the exact complete N strings first found by Baxteiil where 

v k = v Q -2ihy for k = 0,1,---,N- 1. (1.35) 

In ref.i we studied these string solutions of Bethe's equation ( |1.21| ) in the limit as 7 — > 70 
and found a set of equations which determines the parameter vq of the exact complete N 
string ( [L.35| ) in this limiting sense. However, this is not the only way to proceed because 
in the direct coordinate space solution for the eigenfunctions it was shown by Baxter inH 
that there are certain equations which vanish automatically independent of vq and hence 
solutions of the string form (|1.35|) exist where Vq is arbitrary. It is thus to be expected that 
this arbitrary parameter should be closely related to the variable z in the current operator 
( |1.20j ) of the affine Lie algebra. 

This phenomenon of the automatic vanishing of certain constraint equations in the co- 
ordinate space Bethe's ansatz is mirrored in the algebraic Bethe's ansatz by the vanishing 
of the operatoiS 

JV-l 

JJ B(v - 2ik-f ) = (1.36) 
which by (|1.31j) should create the exact complete N string (|1.35|) . 



We denote the operator which replaces the vanishing string creation operator ( 1.36| ) when 
the root of unity condition (|1.4j) holds by B'^(u), It has the property 

In In 

t(v) n n >= (-i)^ m) n ^wto n > ( l37 ) 

j=l k=l j=l k=l 

for all v,v'j and integer I > 1. We find ( |L37| ) is satisfied for B^ N \v) of the form 

B (N \v) = J: (l[B(v - 2U l0 )) (b,(v - 2zfc 7 o) + X{V ~^ o) B v {v - 2iky )) 
k=0 \l=0 /V Y \ v ) / 

AT-1 \ 

IJ B{v-2ih )\ (1.38) 
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where B y (v) and B v {v) specify derivatives of B{v) with respect to 7 and v respectively. This 
operator, with X(v) and Y(v), arbitrary satisfies the commutation relations 

[B {N) (v),B(v')}=0 (1.39) 
[5W(«),BW(i/)]=0 (1.40) 

for all values of v and v '. 

The functions X(v) and Y(v) depend on the highest weight vector on which the operator 
acts. Denoting by v k the solution of Bethe's equation (|1.21|) which specifies the highest 
weight vector we find 

X ( v )=2iY j /sinh L >-(2/+l)* 7 o) 

^ nLi sinn \{v -v k - 2il^ Q ) sinh \{v - v k - 2i(l + 1)70) 

and 

Ytv)=T sinh "<"-' 2 ' + 1 »^) (142) 

3 IK.1 sinh K« - »* - 2<i Tl>) sinh \( v ~ v " ~ 2 'C + l ho) ' 

Because the Chevalley generators S ±( - N ^ and (anti) commute with T(v) the mode 

operators e(n) and f(n) and their generating functions will (anti) commute with T(v) and 
will therefore obey ( |Q7| ). We propose that the solution B^ N \v) ( |L38|) of (|L37l ), ( |Q9|) 
and ( |1.40|) is proportional to the current E~(z). Therefore the zeroes of the function Y(v) 
are the roots of the Drinfeld polynomial computed in a fashion totally independent of the 
definition ( |1.17|) . 

We will prove in section two that B^ N \v) satisfies ( |1.37| ). In section 3 where we will 
study the polynomial Y(v) and compute its degree for the various classes of states specified 
by 

S z = r (mod N) (1.43) 

studied in refs.0-i. We will also see that this polynomial leads to sum rules for each value of 
S z which are not only valid for roots of unity ( |1.4| ) but which hold for all real values of 7. 
Examples of these sum rules are 



S z = : Y. v k = 0) or ™ ( L44 ) 
fc=i 

■j— 1 ■§■— 1 \— 1 — — 1 

1 : {L — (e 4 ' 7 + e" 47 ) £ e Ufe } JJ e"* = {L - (e 4 ' 7 + e" 4 ' 7 ) £ e"**} [] e - "* (1.45) 

fc=i fc=i 

S z =2: L(L - 1) - 2L(e 47 + e" 47 ) £ 



fe=l fc=l fc=l k=l 



k=l 

4-2 #-2 



+ (e 2i7 + e" 2i7 ) X! e 2l,fc + (e 47 + e" 47 ) 2 (]T e^) 2 J] e " fc 

fc=l A-l / k=l 
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L(L - 1) - 2L(e i7 + e" i7 ) £ e" 



V 



+ (e 2i7 + e" 2i7 ) J]) e" 
fe=i 



fc=l 



f-2 



■2^ + ( e n + e-i-yfi 



fc=l 



n 

fc=l 



;i.46) 



The sum rule ( [L.44 ) is the XXZ version of a sum rule in the XYZ model^ and has recently 
been studied numericallyilj!. The sum rules for S z ^ have not previously been seen. 

In section 4 we discuss the relation between these algebraic Bethe's ansatz computations 
and the theory of finite dimensional representations of affine Lie algebras and quantum 
groups at roots of unity. Concluding remarks are in sec. 5. 



II. DERIVATION OF THE OPERATOR B^ N \V) 

We begin the study of the operator B^ N \v) given by ( |1.38|) by proving the commutation 
relations ( |1.39| ) and ( |1.40| ) which do not depend on the form of X(v) and Y{v). We will 



then use the defining property (|1.37|) to determine the functions X{v) and Y(v) given by 

( CT ) and org). 



A. Proof of ( |TT39| ) 

Write the operator B^ N \v) of ( |1.38| ) as 

= B[ N \v) + B { 2 N \v)/Y(v) 

with 



N /fc-1 \ / N-l 

b[ N \v)= Ml B ( v - 2l ho) (B 7 (v - 2ik lo )) W B(v- 2il lo ) 

k \l=0 ) \Z=fc+l 



B™(v) 



N /fc-1 \ / N-l 

]T \J{B{v-2U lo ) \{X{v-2ik lQ )B v (v-2ik lQ ))[ B{v-2il lo 

k \Z=0 / \^=* : + 1 



Then by differentiating ( |1.25|) with respect to v we find 

[B v (v),B(v')] = 

and hence 

[Bi N \v),B(v')}=0. 
It also follows from (|1.25| ) that 



JV-l 

li[B(u l ),B(v')] = 0. 

1=0 



(2.1) 

(2.2) 
• (2-3) 

(2.4) 
(2.5) 

(2.6) 
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Differentiating this with respect to 7 and then setting ui = v — 22Z70 we find 



[B[ N \v),B(v')] + [Y[ B(v-2ti l0 ),B,(v')}=0. 



N-l 



(2.7) 



1=0 



Both terms in the second commutator are zero because of the vanishing condition (|1.36 ) 
and thus 



[B[ N \v),B(v')]=0. 



(2.8) 



Combining (|2.5|) and (2.8) we obtain the desired result ( |1.39|) . 



B. Proof of flOpp 



Using the decomposition (|2.1|) we see that the desired commutation relation Q1.40T ) will 
follow from the three separate commutation relations 

[Bf\v),Bf\v')} = for j,j' = 1,2 

To prove (2.9) for j = f = 1 we note that from (|1.25|) it follows that 



(2.9) 



N-l N-l 

h[b( UI ), n%)i=o 

1=0 1=0 



and thus 



j2 iV-l N-l 

' r n%),n%)]=o 

1=0 1=0 



di 2 



(2.10) 



(2.11) 



Then if we first carry out the differentiations, then set u\ — v — 2iZ7o, Wi = w — 2iljo and 



then use the vanishing condition ( |1.36|) we see that the only non-vanishing terms are those 
where one derivative acts on a B{u{) and one acts on a B(wi) and thus it follows that 

[fW,fWH0, (2.12) 
To prove ( |2.9|) for j = 1, j' = 2 we differentiate (|2.10|) with respect to 7 and wi to find 

j2 N-l N-l 

= pu ( n b (™i) ) [ n B(vn) ] = 



(2.13) 



w=o 



v l=k+l 



Then as before carry out the differentiation with respect to 7. set v \ = v — 2ilj , wi 
w — 2i/7o, and use the vanishing condition ( |1.36| ) to obtain 



/fc-l \ / N-l 

[b[ n \v), n%) )B w (w k )[ n %)|i=o 

\l=0 ) \Z=fc+l 



(2.14) 



from which ( |2.9| ) follows for j = 1, j' = 2. 

Finally we note that fl2.9|) for j = j' = 2 follows immediately from ([L.25 ). Thus we have 
proven the desired commutation relation ( |1.40| ). 
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C. Derivation of X(v) and Y(v) 



The functions X(v) and Y(v) are uniquely determined by the defining property of the 
operator B^ N \v) ( |1.37|) with I = 1. We begin the derivation by proving the following four 
relations: 

Proposition 



JV-1 /JV-1 

(1) A(v)B v (v') J] B(v' - 2il lo ) = {-l) N+m B v {v') J] B{v' 
1=1 \ 1=0 



2ih ) A(v) 



■1) 



N+m 



2 sinh — v) 



B(v) 



i=i 



1) 



N+m 



2 sinh -(?/ — v — 2i7 ) 



/ JV-1 

IJ - 2z7 7 „) | 

V 



IJ B(v' - 2il l0 ) A(v' - 2i l0 ) 



(2.15) 



5*1 



JV-1 /JV-1 

(2) D(v)B v (v') I] - 2z/ 7o ) = (-l)^™^') J] B W 

1=1 \ 1=0 



- 2ih ) D{v) 



-1) 



N+m 



2 sinh ~(y' — u) 



(n-i \ 
B(v) HB(v'-2iho)\D(v') 

\ l ;° J 



(-i) 



N+m N-l 



\ 



2sinh§(v' -v - 2(JV - l)i7 ) 



i=i 

v#JV-l 



- 2(iV - l)i7o) (2.16) 



(3) A(v)b[ n \v') = (-l) N+m B{ N \v')A(v) 

\N+m ( N-l 



N-l 



+ i v 

^ 2smhUv' -v-2ikj ) 



B(v) 



k=0 



Yl B(v' - 2il l0 ) A(v' - 2ifc 7o ) (2.17) 



i=0 



and 



(4) D(v)Bf \v') = {-l) N+m Bf\v')D(v) 

N-l ( -i\N+m I N-l 



V i 

^ 2smh±(v' -v-2ik l0 ) 



B(v) J] B ( v ' - 2il ^o) D(v' ~ 2ik-f ) (2.18) 



!=0 



Proof 



We begin by iterating the basic relations of the algebraic Bethe's ansatz (|1.25|) - (|1.27|) to 
obtain 
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n / n \ n 

A(v) [I = A IJ A (v) + B(v) J2 A fe 

\l=l / k=l 



1=1 



(2.19) 



(=1 



and 



n = a n b ( u i) d ( v ) + b ( v ) e A J n s «> d k) 



(2.20) 



VJ=1 



k=l \ 1=1 



where 



A = n/( u_ ^) ; h k = g(u k -v) ]Jf(u k -Ui) 
i=i i=i 

n n 

A = n/( M ' -W )' A fc = - U & ) JJ - M fe ) 



(2.21) 
(2.22) 



To prove (|2.15 ) we start with ( |2.19| ) and differentiate with respect to U\ to obtain 

n 

A{v)B v (ux) JJ-B(mj) 

= a^( Ui ) f n ^ + ^(^)Ai fn Mm) 



\l=2 



\l=2 



-B(v)J2^kB v ( Ul ) [llBim) A(u k ) 



k=2 



1 = 2 



^((nsw)^) + Bw|;ff[g«wj^) 



(2.23) 



where 



dAi 
9m i 

<9A fc 
9m i 



-/„(v -ui) ]Jf(v -ui) 

1=2 

n n n 

9v{Vl - «) I! ~~ M ') + ~~ W ) J2 fv( U l - U k ) Yl f(ui - Ui) 



1=2 



k=2 



1=2 



-g(u k - v)f v (u k - Mi) J| /(ujfe - vi) for fc ^ 1. 



(2.24) 
(2.25) 

(2.26) 



1=2 



Now after differentiating set n = N and u\ = v' — 2i(l — l) 7 o- Then 

/(-2i 7o )=/(2(iV-l)i T o) = 
/«(-2i7 )=/ t ,(2(JV-l)i 7o ) 



1 



2 sinh i 7o 



(2.27) 
(2.28) 
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and therefore 

iV 



A' 



A=n/(«-«i)=n 



gAi 
du\ 
dA 2 
du\ 
dA k 
du\ 



sinh \ {v — v' + 2iZ7 ) 
; , t Ji sinh |(u - t>' + 2i(Z - 1)70) 

for jfe = 1, • • • , N 



-1 



siiih |(v - «' + 2iJV7 ) 



sinh Ay — v') 



2 sinh -(?/ 



2 sinh !(i>' — i> — 2«7o) 
for k — 3,- ■ ■ .N 



(2.29) 
(2.30) 

(2.31) 
(2.32) 
(2.33) 



and using the vanishing condition ( 1.36|) we find that the only terms in ( |2.23|) which do not 
vanish are the first term and the last term with k = 1,2. Thus we find that (|2.23|) reduces 
to fl2.15|) as desired. 

To prove Q2.17|) we similarly differentiate (|2.19|) with respect to 7 and then set n = N 
and Ui — v' — 2(1 — 1)270- the derivative of the left hand side of ( |2.19|) is 



N /k-1 



N-l 



A (v) E ( II B ( v ' - 2 ^o) ) B,(v' - 2ik lo ) I H B ( v ' - 2 *ho) ] = A(v)B?\v') (2.34) 



k=l \i=o 



A=k+1 



where the term with A 1 (vo) vanishes because of the condition Ql.36 ). To differentiate the 
right hand side of ( [2.19 ) we use both the vanishing condition ( |1.36|) and conditions (|2.29 ) 
and (|2.30|) on A and Afc. Thus we find that the derivative of ( |2.19|) with respect to 7 reduces 
to 



A(v)b[ N \v'^ 



l Br>(v>)A(v) + B(v)j:^ 

k=l 1 



^N+m r>(N)f i\ 



l[B(v'-2iho) A(v'-2ik lQ ). 

(2.35) 



By differentiating ( |2.21| ) with respect to 7 and noting that 



we find 



.cosh ~(v + 2ry) 

% 

sinh(t>/2) 

/ |\JV+m+l 



(2.36) 



$7 sinh \(v' — v — 2i(k — 1)70) 



^2.37) 



we find that ( |2.35| ) reduces to ( |2.17| ) as desired. 

The proof of (|2.16|) and ( |2.18|) follows in an similar manner by differentiating ( |2.20| ) 

QED 
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We now use (FlgD -flngD, the definition of BW(«) in terms of B^u) and 
to find 



A{v)B (N \v') 



-l) N+m B( N \v')A(v) 

. B (v) J] ( X ^'- 2 ^ k ~ ^o) ~ X ( v ' ~ 2ik ^) 

k=0 V 



TV 



l[B{v'-2il l0 ) 



1=1 



Y(v') 

(-1)^+^(^-2^70) 
2 sinh — v — 2ik^ Q ) 



2i 



(2.38) 



D{v)B {N \v') 



-l) N+m B^ N \v')D{y) 

B(v) f X(v'-2 i (k + l) l0 )-X(v'-2tk l0 ) + ^ 

k=0 \ 



X 



UB(v'-2li l0 ) 



(=1 



- 2^70) 
2 sinh \ {v' — v — 2ik r yo) 



(2.39) 



These two operators may now be applied to the Bethe vector ( |1.31| ) . Then we use ( |2.19|) and 
( p.2Q ) to move A{v' — 2ik^o) and D{v' — 2ikj ) to the right to act on |0 >, use the vanishing 
condition ( |1.36|) on all terms generated which contain B(v' — 2i/c7o) and then add the two 
expressions together using the definition of the transfer matrix in terms of A(v) and D(v) 
( |Q4p to find 

n n 

T(v)BW(v') J] B(v k )\0 >= (-l) N+m B^(v')T(v) J[ B(v k )\0 > 

k=l k=l 

~B(v) £ I ( ^ ~ 2%{k ~ 1)70) ~ X{V ' ~ 2lkl ° ] ^ 



fc=0 



N 



X 



+ 



X 



n B(v'-2il l0 ) 



1=1 



Y(v>) 

(-l) N+m a(v' - 2ikj ) 
2 sinh \ {v' — v — 2ikjo) 



2i 



'X(v' - 2i(k + 1)70) - X(v' - 2ik<y ) 



Y{v') 



+ 2% 



( 



Y[B(v'-2U l0 ) 



1=1 



(-l) N+m d(v' - 2ik l0 ) 
2 sinh \ {v' — v — 2ik^ Q ) 



\{B(v k )\Q> 



(2.40) 



k=l 



where 



a(v)= sinh L -(v - 270) ]J f( v ~ v k) 
1 k=i 



d(v)= sinh L -(v + i<y ) Y[ f( v k 
1 k=i 



v . 



(2.41) 
(2.42) 



In order for the defining relation for the operator B^ N \v) to hold for I = 1 the term 
proportional to B(v) in ( |2.40| ) must vanish. Thus we find that X(v) and Y(v) must satisfy 
the following functional equation for k = 0, • • • , N — 1 
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'X(v - 2i(k - 1)70) - X(v - 2ik^ ) 



Y(v) 



— 2i a(v — 2ik^o) 



+ 



'X(v - 2i{k + 1)70) - X(v - 2ik^ ) 
YM 



+ 2% d(v - 2ik^ ) = 0. 



(2.43) 



It is easily verified that this equation (and the periodicity condition X{v + 2iNj ) = X(v) 
is satisfied if 



X{v) X(v - 2i 7o ) 



Y(v) 



Y(v) 



-2i 



Na(v) 



Elo 1 a(v-2zh Q ) 



with 



a(v) = sinh L -(v - 270) f[ - 



\ sinh \{v k - v) sinh \{v k — v + 2ry c 



^2.44) 



(2.45) 



The recursion relation (|2.44j ) only defines X(v) up to an additive function with the 
periodicity X Q (v — 2i^ ) = X (v). But any such function X (v) gives a vanishing contribution 
when used in the expression for Bw(«) ( |1.38| ) by use of the relation which follows from 
vanishing condition (|1.36|) that 



7V-1 



N-l 



£ Bv(v - 2ik l0 ) [] B(v - 2U l0 ) = 0. 



(2.46) 



k=0 



1=0 



Therefore we easily verify by substitution that the only solution ( |2.44 ) which gives a non- 
vanishing contribution to B^ N \v) is 



Af-l 



X(v) = 2iJ2 ka(v -2ik-f ). 

k=0 



Therefore it also follows that 



N-l 

Y(v) = "£a(v- 2ik lQ ) 

k=0 



(2.47) 



(2.48) 



and thus ( |1.41| ) and ( |1.42[ ) are proven. 

Finally we note that Y(v) has a close connection with the eigenvalues of the transfer 
matrix. To see this we note that from ( 1.34| ) we find 



Y\k=i sinh \{v-v h + 2ry ) sinh \{v - v k - 2ry ) 



sinh L \ {v — ij ) 



+ 



sinh L \(v + z7o) 



. (2.49) 



rifc=i sinh \{v - v k - 2ry ) sinh \(v - v k ) IILi sinh \(v - v k + 2«7 ) sinh \(v - v k 

Thus if we replace v by v — 2U^q and sum on I from to N — 1 we find that the two terms 
on the right hand side are equal by periodicity and thus we obtain 

N-l 



1 



t(v — 2il^ ) 



L=0 



U k =i sinh \{v - v k - 2(1 - l)ij ) sinh \(v - v k - 2(1 + 1)170) 



Y(v) (2.50) 



where we have used the expression for Y(v) (\1A2\) . 
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D. Multiple excitations 



It remains to demonstrate that the operator B^ N \v) which has been shown to satisfy 
the defining equation ( |1.37| ) for I = 1 in fact satisfies ( |1.37| ) for all integer / > 1. Consider 
first I = 2. Then using (|2.40| ) to commute T(y) to the right past B^(v[) we obtain 



T(v)B^ N \v[)B (N \v' 2 ) J] B(v k )\0 >= 

k=l 

n 

+ (-l)^B^\v[)T(v)B^(v' 2 ) J] B(v k )\0 > 

k=l 

'X(v[ - 2i{k - l) 7o ) - X{v[ - 2ik lo ) 



N-l 

-B(v) £ 

k=0 

N 



X 



+ 



n^K-2z/ 7o ) 



1=1 



Y(v[) 

2 sinh h(v[ — v — 2ik^o) 



2i 



BW(v' 2 )l[B(v k )\0> 



k=i 



'IK - 2i(k + 1)70) - X(v[ - 2ikj ) 



Y(v[) 



2% 



\{B(v l l -2U lQ ) 



(-l) N +™D(v[-2ikj ) 
2 sinh \ {y [ — v — 2z/c7 ) 



B^(v' 2 )l[B(v k )\0> 



(2.51) 



k=l 



In the first term on the right hand side of (|2.51| ) we use (|2.40|) once again to (anti) 
commute T(v) past B^ N \v 2 ) without picking up any additional terms. In the remaining 
terms we commute A(v[) and D(v[) to the right past B^ N \v' 2 ) using (|2.38|) and (|2.39| ). The 
vanishing condition ( |1.36| ) shows that only the first terms in ( |2.38| ) and ( |2.39| ) contribute 
and thus we find 



T(v)B^\v[)B^ N \v' 2 )l[B(v k )\0>= 

k=l 

+(-l)^ N+m ^B^(v' 1 )B^ N \v' 2 )T{v) f[ B{v k )\0 > 

k=l 

Wl ( 'X{v[ - 2i{k - l) 7o ) - X(v[ - 2ik l0 ) 



+ (-l) N + m B(v)B {N \v' 2 ) 
( 



k=0 



Y(v[) 



-2% 



x 



+ 



N 



Y[B(v' 1 -2U 7o ) 



^ N+mA ^- 2 ^f[B( Vk) \0> 



2 sinh 2 (v[-v- 2ik^) k=1 
X(v[ - 2i(k + l) 7o ) - X(v[ - 2ik lQ ) 



i=i 



Y(v[) 



+ 2i 



( 



Y[B(v[-2li l0 ) 



(2.52) 



2 sinh 2 (v[-v- 2ik^ ) fc=1 
The terms to the right of B^ (y' 2 ) are now exactly the terms shown to vanish in the previous 
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section and thus ( |1.37| ) holds to I = 2. Iteration of this argument proves (|1.37 ) for arbitrary 
I. 



III. PROPERTIES OF X(V),Y(V) AND THE SUM RULES 

The operator B^ N '(v) has the possibility of having poles at the poles of X(v) or the 
zeroes of Y(v ). 



A. The poles of X(v) 

The function X(v) ( |1 .41 ) has poles at 



v = Vj + 2ip r y Q where p = 0, • • • , iV — 1 (3.1) 
As v — > Vj + 2ip7o for p > 1 X(v) behaves as 



X(v) 



2ipsinh L \(vj — 270) 



sinh \(vj + 2ip7o - v) Ut^j sinh |(«* - v j) life sinh - vj + 2ry ) 

+ 2i(p-l)smh L l( Vj +t l0 ) 

sinh |(u,- + 2zp7 - u) sinh \(v k - Vj) Uk sinn \{ v k ~ Vj ~ 2«7o) * 



which upon using the Bethe's equation ( 1.21 ) reduces to 



X( V ) -> -2zsinh L ifa +^o) (3 3) 

sinh i (uj- + 2i/7 - u) Uk^j sinh 5 - v j) Uk sinh I ( u fe - v j - 2 ^7o) 



As v — > f j we have 



w v) ^ 2z(iV-l)sinh L fa + 27o) ^ (3 4) 

sinh \ (vj - v) Ut^j sinh \ i v k ~ Vj ) Ut sinh 5 - v j ~ 2i lo) ' 



Using these limiting forms in the expression for the operator B^ N \v) and using ( 2.46 ) it 
follows that as v — > Vj + 2il^ for any I that 

B (N)( v j K 2iNsmh L \{y j + ij ) 



sinh \{vj - v) TTfe^j sinh \{v k - Vj) U k sinh \{v k - Vj - 2ry ) 

N-l 

xB v { Vj ) [] B( V j-2U l0 ) (3.5) 
1=1 

This operator vanishes when applied to the Bethe state because of the vanishing condition 
( |1.36j ) and therefore B^ N \v) does not have a pole at v = Vj + 2ipryo even though X(v) does. 
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B. Properties of Y[y) 

The function Y(v) denned by ( |1.42 ) also appears to have poles at v — > vj + 2ip~f . In this 
limit 

Y f v ) ^ sinh L |fa -^ ) 

" ' ' ' sinh \{vj + 2ip7o - v) Uh^j sinn l( v k ~ Vj) lit sinh \{v k - Vj + 2ij ) 

+ sinh L ^fa+i7o) 

sinh \{vj + 2ip7 - v) JJk^j sinh \{v k - v-) Hi sinh \(v k - Vj - 2i^ ) ' 

However by use of the Bethe equation (|1.21|) we see that the residue at this pole vanishes. 
Therefore Y(v) in fact has no poles and therefore must be a Laurent polynomial in e v . 
Moreover it follows from the periodicity condition Y(v + 2270) = Y(v) that Y(v) is in fact 
a Laurent polynomial in 

z = e Nv . (3.7) 

To further study the properties of the Laurent polynomial we introduce the positive and 
negative degrees of the polynomial defined by the behavior of Y(v) as v — > ±00 as 

Y{v) ~ C±e ±Nd±v as^ ±00. (3.8) 

We may thus define the Drinfeld polynomial as 

Pa(z) =e d - Nv Y(v). (3.9) 

This is a polynomial in z with degree 

d = d + + d_ (3.10) 

This limiting behavior of Y(v) is obtained from ( |1.42|) . There still are several cases to 
consider depending on the number of Bethe roots v k which take on the value of ±00. 
Consider first the case of no infinite v k . Then as v — > ±00 we find that 

AT-l N-l 

Y(v) ~ 2- {L - 2n) C±(l)e^~ n - l) e~^' n - l)1 ^ i2j+1) 
1=0 j=0 

= 2^ L - 2n ^C ± (r)Ne lv ^- n ~ r) (3.11) 



where 



^ - n - r = (mod N) and r = 0, 1, • • • , N - 1 (3.12) 



and the first few C±(j) are 
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c±(o)= n e±vk ( 3 - 13 ) 
fe=i 

(n \ n 

-L + (e i7 ° + e- i7 °) ]T e ±t,fc J] e ±Vk (3.14) 
fe=i / fe=i 

C±(2)= ( - 1) - L(e* 70 + e"* 70 ) £ e ±l "= 

V 2 fc=i 

-r n i n \ n 

+ i( e 2no + e -2i 70 ) £ e ±2v k + V e i 70 + e -^) 2 (^ e ±,;fc ) 2 J] e ±,;fc (3.15) 
fc=i fc=i / fc=i 

Therefore for the case of no infinite roots Vk we have 

d + = d„ = (^-n-r)/N=[S z /N} (3.16) 

where in the last line we have used ( |1.33|) and [x] denotes the greatest integer in x Thus the 
degree of the Drinfeld polynomial for highest weight states with no infinite roots is 

d = 2[S z /N}. (3.17) 

When S z = (mod N) all evidence in refs.0-i is that infinite roots never occur. However 
for all other cases there are indeed infinite roots and the evidence of refs.i-i is that there are 
two cases to distinguish 

(A) The pair v = oo and — oo occurs p times. 

(B) The single root v k = oo or — oo occurs s times. 

For case (A) we let p pairs of roots go to ±oo in ( |1.41| ) and ( |1.42| ) and define 

Y(v) Y p {v) 

where X p {y) and Y p (v) are obtained from ( |1.41|) and ( p..42| ) by merely omitting the infinite 
roots (and thus replacing n by n — 2p). 

The previous argument generalizes by replacing n by n — 2p and thus we find 

Y p (v) ~ 2 -^ 2n+4 ^C ± (r)Ne^~ n - r+2p) (3.19) 

where 

^—n — r + 2p = (mod N). (3.20) 

From the data of table 10 of ref.i for iV = 3 we find that for S z = 1 (mod 3) that p = 2 and 
for S z = 2 (mod 3) that p — 1. This gives r = p for iV = 3. If this holds generally for all iV 
then we find from ( |3.20| ) that the degrees of Y p (v) are 

d + = d_ = [S*/N] + 1 (3.21) 

and thus the degree of the Drinfeld polynomial is 
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d = 2[S Z /N] + 2 (3.22) 
For case (B) we let s roots Vk — > oo or — oo in ( |1.41|) and ( 1.42 ) to find 



and 



/e^( 2 ' +1 ^» sinh L - (2/ + l)^ 7o ) 
(v) ±s - 2z ^ n n-. sinh i (v _ Vk _ 2Ul0 ) sinh I( W - v k - 2i{l + l) 7o ) (3 " 23) 



^ e^ S inh^(.- i7o (21 + l)) 

to IK sinh |( Ufc - u + 2il lo ) sinh ±( Ufc - v + 2i(l + l) 7o ) 1 ' ] 



Now the limiting behaviors as v — > ±oo are no longer the same and we find from ( |3.24 ) that 
as f — > oo 



Y± s (v) ~ 2-( L - 2n - 2s )C±(r + )iVe |,;|( *- n - r++s) (3.25) 



where 



^-n-r + + s±s = (mod JV). (3.26) 

and similarly as u — > — oo 

y± a (v) ~ 2 - {L - 2n - 2s) C±(r_)Ne lv ^~ n - r - +s) (3.27) 

with 

^-n-r_ + s=Fs = (mod N). (3.28) 

In ref.i we found that for S z = 1 (mod 3) that s = 2 and if 5* z = 2 (mod 3) that s = 1. 
From this we find for Y +s if S z = 1 (mod 3) that r + = 2, r_ = 1 and if S* 2 = 2 (mod 3) that 
r + = 1, r_ = 2. Thus we find for this (and for all other N with the same relation of r± to 
S z that 

d + = [S z /N] + 1, d_ = [S z /N] (3.29) 
and thus the degree of the Drinfeld polynomial is 

d = 2[S z /N] + l. (3.30) 



C. Sum rules for Bethe roots 

In general the Laurent polynomial Y(v) has d+ + gL > 1. However if the highest weight 
state has no infinite roots and if 

0<S Z <N-1 (3.31) 
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then we see from ( |3.11| ) that Y[v) is a constant which must equal the limiting values obtained 



from v — > ±00. Thus it follows that when (|3.31|) holds we have N distinct sum rules 

C + (S Z ) = C^{S Z ). (3.32) 
For example we find from ( p,13[ )-( p.l5[ ) that 



L/2 

S z = : Vf~ = 0, or Tii (3.33) 



k=0 



-1 



S z = 1 : {L — (e i7 ° + e~* 70 ) £ e"*} ]J e Vk = {L - (e^° + e" i7 °) £ e^*} f[ e""* (3.34) 



=1 fc=l k=l k=l 



1 



S z = 2 : -L(L - 1) - L(e i7 ° + e" i7 °) ^ 

V 2 fc=i 



1 n 1 n \ n 

+ _( e 2no + e -2i 7 o) ^ e 2v k + (e no + e- i70 ) 2 (^ e" k f J] e " fc 
fc=i 2 fe=i / fc=i 




+ I( e 2no + e -2i 7 o) ^ e -2 Vfe + i ^ e i 7 o + e- 470 ) 2 (^ e^ fc ) 2 ) J] e~ Vk (3.35) 
2 fc=i 2 fc=i / fc=i 

These sum rules have been derived under the assumption that 70 satisfies the root of 
unity (rationality) condition ( |1.4| ). However, for fixed S z , these sum rules hold for all N such 
that S z < N — 1 and will therefore hold for all 7 such that 7/71 is irrational. Moreover the 
sum rule will also hold at all rational values of j/n for the Bethe's roots which are obtained 
by continuity from the irrational values of j/n provided that all the Vk are finite. Thus we 
have proven the general sum rules ( |1.44| )-( |I.46| ) given in the introduction. 



It is worth noting that these sum rules follow merely from the properties of the function 
Y(v) as defined in terms of the Bethe roots by ( |1.42j ). The relation of Y(v) to the operator 
B( N >(v) has not been used in this derivation and in principle these sum rules, which to the 
authors knowledge are new, could have been derived without the knowledge of the sl 2 loop 
algebra symmetry which forms the starting point for this present paper. 

IV. REPRESENTATION THEORY 

The mathematical theory of finite dimensional representations of affine Lie algebrasS-@ 
and the demonstrationlll of the sl 2 loop algebra symmetry in the six vertex model at roots 
of unity in the sector S z = (mod N) makes no use of the algebraic Bethe's ansatz. On 
the other hand the computation of the B^ N '(v) in this paper makes no use of representation 
theory. Nevertheless the conjecture that 

B {N \v) = N(z)E-(z), (4.1) 

where N(z) is a scalar with no poles implies that the two methods are part of the same 
subject and that by combining them we are able to obtain results which have so far been 
inaccessible to either method in isolation. 
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One such result is the explicit expression for the Drinfeld polynomial 

Pa{z) =e d - Nv Y(v). (4.2) 

Furthermore if the roots of Y(v) are distinct the current E~(z) will be given solely in 
terms of the residue at the poles z = a,j where Pn(aj) = 0. Therefore is must follow that 

When the roots are distinct representation theory says that only spin 1/2 representations 
occur in the degenerate multiplets and therefore 

E n N)2 (j)\n>=0. (4.4) 

From this we see that if we denote the spin in the multiplet by 

S z = S^ ax - Nl with < / < d n (4.5) 

that the multiplicity is given by the binomial coefficient given in the introduction ( |1.19|) 



multiplicity = ( ^ J . (4.6) 

This is in complete agreement with ref .0 for all sectors not only for S z = (mod N) where 
the symmetry algebra was proven but in all other sectors where a projection was needed in 
order to obtain the algebra. Thus the mechanism of the infinite roots which automatically 
appears in the algebraic Bethe ansatz makes explicit the projection mechanism only effected 
on the computer when using the Chevalley basis for the loop algebra generators^. We also 
note from (|4.3| ) that the operator B^ N '(v) acting on \Q > will create a vector space of 
dimension dn. For the case N = 3 this dimension has been computed from the coordinate 
space Bethe's ansatz by Braak and Andrefl 



V. CONCLUSIONS 

The results of the previous section are obtained by using the representation theory of 
affine Lie algebras§-0 to provide existence theorems for algebraic Bethe's ansatz compu- 
tations. It is clearly desirable to avoid representation theory altogether and to explicitly 
compute the normalizing constant N(z) and the residue operators E n N \j) explicitly from 
B( N \v). Conversely it is desirable to compute the expression (|1.42|) for the Drinfeld poly- 



nomial directly from representation theory without the use of the algebraic Bethe's ansatz. 
These can perhaps be viewed as mathematical problems. 

But from the physical point of view perhaps the most interesting question is to inquire 
into the physical meaning and significance of the Drinfeld polynomial ( 1.42Q . This polyno- 



mial has been seen previously in our study of the limiting form of Bethe's ansatz as 7 — > 70 
(see the function K(aj) of eqn. 2.18 of ref.!). In addition very closely related expressions 
have been seen in the computation of the free energy of the superintegrable chiral Potts 
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model (see eqn. 4.5 of ref.EH) and in the study of the RSOS models (see eqn. 3.29 of ref.Hl). 
The evaluation parameters, which are the roots of the Drinfeld polynomial, must have a 
further physical significance which is yet to be discovered. 

Finally it needs to be noted that everything in this paper can be extended from the spin 
1/2 XXZ model to any model based on a quantum group just as the loop symmetry was 
extended in ref.il. 



22 



Acknowledgments 

We are very pleased to acknowledge fruitful discussions with R.J. Baxter, E. Date, T. 
Deguchi, M. Jimbo, M Kashiwara, C. Korff, A. Kuniba, T. Miwa, and K. Okado. We also 
thank M. Kashiwara and T. Miwa for their invitation to "MathPhys Odyssey 2001 Integrable 
Models and Beyond" where this work began. One of us (BMM) is pleased to thank M. 
Kashiwara and the Research Institute of Mathematical Sciences of Kyoto University for 
their hospitality. This work was supported in part by the National Science Foundation 
under Grant DMR-0073058. 



23 



REFERENCES 



1 T. Deguchi, K. Fabricius and B.M. McCoy, The sl^ loop algebra symmetry of the six 
vertex model at roots of unity, J. Stat. Phys. 102 (2001) 701. 

2 K. Fabricius and B.M. McCoy, Bethe's equation is incomplete for the XXZ model at roots 
of unity, J. Stat. Phys. 103 (2001) 647 

3 K. Fabricius and B.M. McCoy, Completing Bethe's equation at roots of unity, J. Stat. 
Phys. 104 (2001) 575. 

4 V. Chari, Integrable representations of affine Lie-algebras, Inv. Math. 85 (1986) 317. 

5 V. Chari and A. Pressley, Quantum affine algebras, Comm. Math. Phys. 142 (1991) 142. 

6 V. Chari and A. Pressley, Quantum affine algebras and their representations, in Repre- 
sentations of groups, Ed. B.N. Allison and G.H. Cliff, Canadian Mathematical Socisty 
Conference Proceedings 16 (1995) 59. 

7 V. Chari and A. Pressley, Quantum affine algebras at roots of unity, Representation Theory 
1 (1997) 280, |q-alg/ 969003 1| 

8 V.G. Kac, Infinite dimensional Lie algebras 3rd ed. (Cambridge Univ. Press 1990). 
9 H.A. Bethe, Zur Theorie der Metalle, Z. Physik 71 (1931) 205. 

10 R. Orbach, Linear antiferromagnetic chain with anisotropic coupling, Phys. Rev. 112 
(1958) 309. 

11 L.R. Walker, Antiferromagnetic linear chain, Phys. Rev. 116 (1959) 1089. 

12 J. des Cloizeaux and M Gaudin, Anisotropic linear magnetic chain, J. Math. Phys. 7 
(1966) 1384. 

13 C.N. Yang and CP. Yang, One-dimensional chain of anisotropic spin-spin interactions 

I. Proof of Bethe's hypothesis for ground state in a finite system, Phys. Rev. 150, 321 

(1966) ; 

14 C.N. Yang and CP. Yang, One-dimensional chain of anisotropic spin-spin interactions. 

II. Properties of the ground-state energy per lattice site for an infinite system, Phys. Rev. 
150, 327 (1966). 

15 E.H. Lieb, Exact solution of the problem of the entropy of two-dimensional ice, Phys. 
Rev. Letts. 18, 692 (1967). 

16 E.H. Lieb, Exact solution of the F model of an antiferroelectric, Phys. Rev. Lett. 18, 1046 

(1967) . 

17 E.H. Lieb, Exact solution of the two-dimensional Slater KDP model of a ferroelectric, 
Phys. Rev. Lett. 19, 108 (1967). 

18 E.H. Lieb, Residual entropy of square ice, Phys. Rev. 162, 162 (1967). 

19 L.A. Takhtadzhan and L.D. Faddeev, The quantum method of the inverse problem and 
the Heisenberg XYZ model, Russian Math. Surveys 34 (1979) 11. 

20 V.E. Korepin, N.M. Bogoliubov, and A.G. Izergin, Quantum Inverse Scattering Method 
and Correlation Functions, (Cambridge University Press 1993). 

21 R.P.Langlands and Y. Saint-Aubin, Algebro-geometric aspects of the Bethe equations, in 
Proceedings of the Gursey Memorial Conference I-Strings and Symmetries, ed. G. Atkas, 
C. Saclioglu and M. Serdaroglu, Lecture Notes in Physics, vol. 447, Springer Verlag (1994) 

22 V. Tarasov and A. Varchenko, Bases of Bethe vectors and difference equations with regular 
singular points, International Mathematics Research Notes 13 (1995) 637 q-alg/9504011 



2>i 



R.P. Langlands and Y. Saint-Aubin, Aspects combinatoires des equations de Bethe, in Ad- 



24 



vances in Mathematical Sciences: CRM's 25 five years, Ed. Luc Vinet, CRM Proceedings 
and Lecture Notes vol. 11 Am. Math. Soc. (1997) 231. 

24 R.J. Baxter, Eight-vertex model in lattice statistics and one-dimensional anisotropic 
Heisenberg chain III: Eigenvectors of the transfer matrix and the Hamiltonian, Ann. Phys. 
76 (1973)48. 

25 R.J. Baxter, Eight-vertex model in lattice statistics and one-dimensional anisotropic 
Heisenberg chain I: Some fundamental eigenvectors, Ann. Phys. 76 (1973) 1. 

26 V. Tarasov, Cyclic monodromy matrices for sl(n) trigonometric R- matrices, Comm. Math. 
Phys. 158 (1993) 459. 

27 R.J. Baxter, Partition function of the eight-vertex lattice model, Ann. Phys. 70 (1972) 
193. 

28 A. Doikou and R.I. Nepomechie, Discrete symmetries and the S matrix of the XXZ chain, 
J. Phys. A31 (1998) L621. 

29 A. Doikou and R.I. Nepomechie, Parity and charge conjugation symmetries and the S 
matrix of the XXZ chain, in Statistical Physics on the Eve of the 21st Century, ed. M.T. 
Batchelor and Luc T. White (World Scientific 1999) 391. 

30 D. Braak and N. Andrei, On the spectrum of the XXZ-chain at roots of unity, |cond- 
mat/0l06593|. 



31 G. Albertini, B.M. McCoy, and J.H.H. Perk, Eigenvalue spectrum of the superintegrable 
chiral Potts model, Advanced Studies in Pure Mathematics, vol. 19 Ed. M. Jimbo, T. 
Miwa and A. Tsuchiya, (Academic Press, 1989) 1. 

32 V. V. Bazhanov and N. Yu. Reshetikhin, Critical RSOS models and conformal field theory, 
Int. J. Mod. Phys. 4 (1989) 115. 

33 C. Korff and B.M. McCoy, Loop symmetry of integrable vertex models at roots of unity, 
hepth/0104120. 



25 



